Collapse of the critical state in superconducting niobium 
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Giant abrupt changes in the magnetic flux distribution in niobium foils were studied by using magneto- 
optical visualization, thermal and magnetic measurements. Uniform flux jumps and sometimes almost 
total catastrophic collapse of the critical state are reported. Results are discussed in terms of thermomag- 
netic instability mechanism with different development scenarios. 



Large jumps in magnetization, temperature, ultrasonic 
attenuation and resistivity have been observed in many su- 
perconducting materials since early 60s 00 000 01- 
Although no definitive mechanism has been identified, it is 
believed that such changes are caused by sudden redistribu- 
tion of Abrikosov vortices triggered and supported by ther- 
momagnetic instabilities fl UE EE EH El ■ Di- 
rect observations of flux jumps using miniature Hall probes 
13} 1 141 Il5l 1 1611 and real-time magneto-optical imaging 
3L ll7lll8Lll9Ll20L I2H1 have confirmed precipitous propa- 
gation of flux into the specimen. Although some evidence 
of large - scale jumps exists (3, 1 1 3fl _ direct observations 
of avalanche - like behavior in films have mostly found 
dendritic and branched finger-like patterns (on a macro- 
scopic scale of the whole sample) or small jumps (few 
tens of micrometers) forming feather-shaped flux fronts 

HEiElElElEfllll. 

Dendritic avalanches are expected for thin films with 
high critical current density I Vt flllll2l l - The problem is that 
dendritic avalanches are too small to result in the observed 
jumps of magnetization, M. If SM is a variation of mag- 
netization during a flux jump, then the dendritic avalanches 
contribute 6M < 0.1M I19B . However, the observed 
changes are often much larger 1161 12211 (also see Fig. [Q. 
They must be associated with more dramatic catastrophic 
changes of the critical state when a nonuniform (Bean's) 
distribution of vortices suddenly collapses. In this paper 
we present direct observations of the catastrophic collapse 
of the critical state in niobium foils and provide evidence 
for thermomagnetic origin of this effect. 

Pure Nb (99.99%) strips of various thickness (5, 10 and 
25 fim) were obtained from Goodfellow. Sample widths 
varied from 0.7 to 4 mm and the lengths varied from 2 to 
10 mm. Magnetization loops, M (H), were measured in 
a Quantum Design MPMS. Resistivity, specific heat and 
transport critical current were measured in a Quantum De- 
sign PPMS. We also measured sample temperature varia- 
tion in quasi-adiabatic conditions and found temperature 
jumps associated with the flux jumps. 

Figure^shows M (H) in a 25 /im thick Nb foil at dif- 
ferent temperatures, (a) 1.8 K, (b) 3 K and (c) 4 K. At 
the lowest temperature, Fig. ^a), jumps in M (H) are 



observed both upon penetration after cooling in zero-field 
(ZFC) and upon flux exit. At higher temperatures, jumps 
were only observed upon reduction of the magnetic field, 
(b), and no jumps were observed above 4 K, (c). In Fig. 
Qb) and (c) curves are shown after the external field was 
reduced from 5 kOe to zero and an attempt was made to 
observe jumps upon increasing the magnetic field of the 
same orientation. Unlike the case of the opposite field, no 
jumps in M (H) were observed. These results hint at ther- 
momagnetic instability mechanism 00. The comparison 
of flux jumping instability scenarios in thin foils and bulk 
samples a large Nb single crystal (5 mm diameter, 1 mm 
thickness) was measured at several temperatures. Figure^ 
(d) shows M {H) at T = 1.8 K for both ascending and de- 
scending branches. For other temperatures only descend- 
ing branches are plotted. Giant magnetization jumps were 
observed upon decrease of the magnetic field at tempera- 







T = 1.8 K 




T = 3 K 


25 urn strip 


V 


(a) 


25 Lim strip 


(b) 






T = 4 K 


A 1 II 6.5 K > 


5 K 
a/ 5.5 K h 


25 jim strip 




(c) 


' ' ' 1 ' J ,6 K 

1 mm thick \ 
single crystal \ 


/ / 
4.5 K / 

S4*t 



-2000-1000 1000 2000 

H (Oe) 



-2000-1000 1000 2000 
H (Oe) 



FIG. 1: (color online) Magnetization loops measured in 25 /im 
Nb foil at different temperatures, (a) 1.8 K, (b) 3 K, (c) 4 K and, 
for comparison, several M(H) curves measured from 4 kOe to —4 
kOe in large Nb single crystal. For (b) and (c) curves start from 
remanent magnetization, see text for details. 
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FIG. 2: (color online) Penetration of magnetic field in 10 fim 
thick Nb strip after zero-field cooling, (a) 160 (b) 415 (c) 480 Oe. 
Penetration in (b) and (c) was avalanche-like at fields indicated. 
(Real-time video online [25]). 



tures up to 5.5 K. Evidently, variation of magnetization at 
a jump is much larger in a thick sample. 

Visualization of the magnetic induction, B (r), on the 
sample surface was performed in a flow-type 4 He cryostat 
with the sample in vacuum attached to a copper heat ex- 
changer (i.e., cooling conditions are different from mea- 
surements in QD MPMS, Fig0. Bismuth-doped iron- 
garnet heterostructure film with in-plane magnetization 
was used as a magneto-optical indicator |23l|24|l. In all im- 




FIG. 3: Catastrophic collapse of the critical state upon entry of 
negative self-field at T = 3.9 K. (a) 250 Oe (b) 200 Oe (c) 100 
Oe (d) 50 Oe. (Real-time video online |H| ). 



ages, intensity is proportional to the magnitude of B (r). 
Moreover, due to the dispersion of the Faraday rotation for 
different wavelengths, up and down directions of the mag- 
netic induction can be distinguished. In our case, green cor- 
responds to flux out of page and yellow corresponds to the 
opposite orientation (antiflux). This information is useful 
for flux- antiflux annihilation experiments described below. 

In magneto-optical experiments, sudden collapses of 
magnetization were observed only below T ~ 4.5 K for 
all samples. Figure |2] shows penetration of the flux into a 
10 /jm Nb strip at increasing magnetic fields, (a) 160 Oe 
(b) 415 Oe and (c) 480 Oe. The large-scale abrupt jumps 
in Fig. 13 (b) and (c) occurred right before fields at which 
images were acquired (about 405 Oe and 470 Oe, respec- 
tively). In our case of 30 Hz data acquisition rate, the jumps 
seem instantaneous. Real-time dynamics can be viewed 
online ll25ll . 

The largest and most dramatic features occur upon de- 
crease of a magnetic field when self-field antiflux induced 
at the edges due to large demagnetization factor starts en- 
tering the sample. The antiflux annihilates with the trapped 
flux providing additional mechanism to trigger thermal in- 
stability. An avalanche in progress is difficult to stop, be- 
cause the whole sample is in the critical state and an in- 
crease of temperature leads to reduction of the local critical 
current density resulting in further flux redistribution. This 
results in large uniform avalanches when the whole criti- 
cal state literally collapses. In Fig. |3]two such collapses 
are shown, - panels (b) (jumped after (a)), and (d) (jumped 
after (c)). Some intensity banding in Fig. |3] reflects the 
structure of the trapped flux in this sample. 
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FIG. 4: (color online) Profiles of the magnetic flux before and 
after collapse of the critical state upon entry of the negative flux 
generated by self-field. The profiles were measured along dashed 
lines in Fig|5](c) and (d). 



FIG. 5: Entry of the antiflux (from right to left): (a) -320 Oe (b) 
—470 Oe. Note a huge uniform avalanche, (c) —740 Oe. Note 
secondary penetration on top of the flux distribution flattened by 
the avalanche. It indicates that j c returned to its large "cold" 
values. (Real-time video online [ 25]). 



Flux profiles corresponding to the collapse of the criti- 
cal state are shown in Fig. |4] The profiles were measured 
along the dashed lines shown in Fig.[3](c) and (d). Clearly, 
at some magnetic field, the slightest variation results in a 
change of the entire distribution of vortices. The observed 
giant jumps may span more than half the width of the sim- 
ple, which, to the best of our knowledge, has never been 
observed before. Moreover, the collapse washes away both 
the interior critical state as well as any newly penetrated 
antiflux at the edges. When the magnetic field of the oppo- 
site polarity is applied, antiflux annihilates with the trapped 
flux similar to the self -field antiflux discussed above. This 
results in large uniform avalanches as shown in Fig. |5] 
In general, the largest avalanches and sometimes total col- 
lapse of the critical state occurs at small fields upon entry 
of an antiflux. We think that this is because critical current 
density is largest and most strongly field - dependent at low 
fields. 

Although thermal instability seems to be the most natural 
scenario, direct evidence for thermal effects is scarce. Tem- 
perature jumps were detected in Nb-Zr alloys |4]. MgB 2 
films coated with gold showed suppression of the den- 
dritic avalanches 12711 . We report direct measurements of 
the temperature variation associated with flux jumps in Nb 
foils. In this experiment, the sample was suspended by 
four wires on a thin sapphire stage normally used for spe- 
cific heat measurements in a Quantum Design PPMS. The 
arrangement allows for direct reading of the sample ther- 
mometer in almost adiabatic conditions. 

Figure[6]shows field sweeps (all taken from 3 kOe to - 3 
kOe). Temperature jumps were detected at the same range 
of temperatures and fields as magnetization measured in 
the same sample. To emphasize the variation, it was mul- 
tiplied by a factor of 20. Indeed, local temperature, T a , 
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FIG. 6: (color online) Temperature variation upon decrease of 
magnetic field from 3 kOe to —3 kOe in a 25 fim thick Nb foil. 
The temperature variation was multiplied by a factor of 20 for 
clear visual appearance. 



within the "volume" of the avalanche, V a ~ 5 x 10 4 
mm , reaches very large values of order of Then 
the average temeperature measured by the calorimeter is 
given by 5T = T a V a /V t , where total volume of the sam- 
ple, V t ~ 5 x 10~ 2 mm 3 . This gives 5T ~ 0.01T C K, in 
agreement with the direct measurements, Fig|5J 

In general, there are two channels to dissipate the en- 
ergy of the critical state - dynamic and adiabatic 00 01- 
The dynamical channel corresponds to the energy removal 
to a coolant and it is characterized by the dimensionless 
parameter , a = pj 2 c d/h (T c - T ) fl IH 
Here j c is the critical current, h ~ 5 W/cm 2 K, is the 
heat transfer coefficient, p ~ 0.5 pil-cm is the nor- 
mal state resistivity, d - thickness and T - temperature 
of the coolant. The adiabatic channel corresponds to the 
energy absorption by the material itself which is accom- 
panied by an increase of the sample temperature. The 
threshold parameter for thin samples is given as (3 = 
Mwj 2 / [c 2 C v p (T ) (T c - T )] Efl], where w is the sam- 
ple halfwidth. Using the results of direct measurements 
on our foils, j c (T) = 1.5 x 10 6 (1 - T/T c ) A/cm 2 and 
C p (T) = 0.25 (T/T c f J/mol-K, we obtain a ~ 8 and 
(5 ~ 15. Large values of a and (3 indicate that both chan- 
nels can take away only insignificant fraction of the Joule 
heat, which makes flux jumps inevitable, in full agreement 
with our experiments. 

A flux jump in a thin sample develops into either den- 
dritic or uniform pattern depending on the ratio of ther- 
mal and magnetic diffusivities, r = K/4npC, where k 
is the thermal conductivity. lfl2ll Qualitatively, if r is very 
small, the lateral heat diffusion is slow and flux propa- 
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FIG. 7: (color online) Profiles of magnetic induction in 2.88 mm 
wide and 10 /im thick strip upon flux penetration after ZFC. Two 
inner profiles correspond to large uniform jumps shown in Fig. 



gates along narrow channels forming a dendritic pattern. 
For large r the heat diffusion establishes essentially uni- 
form temperature distribution and a flux jump spreads out 
uniformly in all directions. The threshold value is given 
(for a 1) as r c ~ w/8dnP 1 12], where n is the ex- 
ponent in the current- voltage law, E oc j n . For n = 30 
and « = 10 W/Km, we find r ~ 1, while r c ~ 0.03 in- 
dicating that the flux jumps should develop uniformly - as 
we indeed observe experimentally. To quantitatively val- 
idate these estimates we considered development of flux 
jumps by numerically solving the thermal diffusion and 
Maxwell equations, C p dT/dt = jE — h(T — T ) jd and 
dE/dx = —dH/dt. The sample was assumed to be a thin 
infinite strip, where all the quantities depend only on the 
coordinate x across the strip. 

Measured (symbols, corresponding to Fig. and cal- 
culated profiles (solid lines) of the magnetic induction are 
shown in Fig. Q If we om it the Joule term jE, the calcu- 
lations produce the well-known Bean-model profiles 
that perfectly fit the experimental profiles before a jump, 
see the 300 Oe curves. If we start from such a Bean pro- 
file and, keeping applied field constant, introduce an in- 
finitesimal perturbation in T, a flux jump occurs. The curve 
for 415 Oe shows the final flux profile after such a jump. 
Apparently, the agreement with experiment is again very 
good, which further supports the thermomagnetic mecha- 
nism of the observed instability. 

Furthermore, we tried to control the rate of heat conduc- 
tance by introducing a thin mylar film between the heat ex- 
changer and the sample. Both, feather-shaped protrusions 
and large uniform avalanches were observed in this case. 
Moreover, uniform jumps appeared even in ZFC samples 



where they were absent before. These experiments demon- 
strate that reducing the heat link to the thermal bath makes 
the avalanches more likely to happen. 

In conclusion, we observed large scale giant uniform 
vortex avalanches and even collapse of the critical state in 
niobium foils. Direct measurements of temperature vari- 
ation as well as analysis of the magnetic and thermal en- 
ergy balance support thermomagnetic instability scenario. 
Comparable values of the magnetic and thermal diffusivi- 
ties (r ~ 1) result in spatially uniform giant flux jumps. 

We thank Alex Gurevich, Igor Aranson and Yuri 
Galperin for useful discussions and Thomas A. Girard for 
providing the samples and for discussions. Ames Labora- 
tory is operated for US DOE by the Iowa State University 
under Contract No. W-7405-Eng-82. R. P. acknowledges 
support from NSF Grant DMR-05-53285 and from the Al- 
fred R Sloan Research Foundation. 



* Electronic address: prozorov@ameslab.gov 
[1] R. G. Mints and A. L. Rakhmanov, Rev. Mod. Phys. 53, 551 
(1981). 

[2] A. V. Gurevich, R. G. Mints and A. I. Rahmanov, "The 
Physics of Composite Superconductors" ', (Begell House 
Inc., New York, 1997). 

[3] E. Altshuler and T. H. Johansen, Rev. Mod. Phys. 76, 471 
(2004). 

[4] L. T. Claiborne and N. G. Einspruch, J. App. Phys. 37, 925 
(1966). 

[5] B. B. Goodman and M. Wertheimer, Phys. Lett. 18, 236 
(1965). 

[6] M. Levy et al, Phys. Rev. B 2, 2804 (1970). 

[7] I. S. Aranson et al, Phys. Rev. Lett. 94, 037002 (2005). 

[8] B. Biehler et al, Phys. Rev. B 72, 024532 (2005) 

[9] A. L. Rakhmanov et al, Phys. Rev. B 70, 224502 (2004). 
[10] D. V. Shantsev et al, Phys. Rev. B 72, 024541 (2005). 
[11] Y. E. Kuzovlev, cond-mat/0607143 (2006). 
[12] D. V. Denisov et al, Phys. Rev. B 73, 014512 (2006). 
[13] S. S. James et al, Physica C 332, 445 (2000). 
[14] E. Altshuler et al, Phys. Rev. B 70, 140505 (2004). 
[15] P. Esquinazi et al, Phys. Rev. B 60, 12454 (1999). 
[16] V. V. Chabanenko et al, Physica C 369, 82 (2002). 
[17] P. Leiderer et al, Phys. Rev. Lett. 71, 2646 (1993). 
[18] C. A. Duran et al, Phys. Rev. B 52, 75 (1995). 
[19] T. H. Johansen etal, Europhys. Lett. 59, 599 (2002). 
[20] M. S. Welling et al, Physica C 411, 11 (2004). 
[21] R. J. Wijngaarden et al, Europ. Phys. J. B 50, 1 17 (2006). 
[22] D. P. Young etal, Supercond. Sci. Technol. 18, 776 (2005). 
[23] A. A. Polyanskii, D. M. Feldmann and D. C. Larbalestier, in 
"The Handbook on Superconducting Materials" , D. Card- 
well and D. Ginley, editors, UK, pp. 1551 (2003). 
[24] C. Jooss et al, Rep.Prog. Phys. 65, 651-788 (2002). 
[25] http://cmp.ameslab.gov/supermaglab/video/Nb.html 
[26] E. H. Brandt and M. Indenbom, Phys. Rev. B 48, 12893 
(1993). 

[27] Eun-Mi Choi etal, Appl. Phys. Lett. 87, 152501 (2005). 



